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ABSTRACT 

Using a one-zone model for the evolution of dust in spiral galaxies and 
applying the instantaneous recycling approximation to the model equations, we 
investigate dust-to-gas ratio of spiral galaxies. Four processes are considered; 
dust formation from heavy elements ejected by stellar mass loss, dust destruction 
in supernova remnants, dust destruction in star-forming regions, and accretion 
of heavy elements onto preexisting dust grains. The equations describe a 
simplified relation between dust-to-gas ratio and metallicity. The relation is 
independent of star formation rate. By comparing the theoretical relation with 
the observational data of nearby spiral galaxies, we show (i) that the accretion 
process onto the preexisting dust particles must be taken into account for the 
spiral galaxies; (ii) that the efficiency of dust production from heavy elements 
ejected by stars can be constrained by the spiral galaxies with low metallicity; 
(iii) that the Salpeter and Scalo initial mass functions are both consistent with 
the data of the spiral galaxies. 

Subject headings: galaxies: evolution — galaxies: ISM — galaxies: spiral - 
ISM: dust, extinction 

1. INTRODUCTION 

Interstellar dust is composed of heavy elements made and ejected by stars. Supernovae 
(SNe) govern the dust content in galaxies: Observationally, dust has been clearly detected 
in the SN1987a event (Moseley et al. 1989; see also Kozasa, Hasegawa, & Nomoto 1989), 
and theoretically, Dwek & Scalo (1980) demonstrated that SNe are the dominant sources for 
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the formation of dust grains. Dwek & Scalo (1980) also showed that the dust is destroyed 
in SN remnants (SNRs). 

Recently, galaxy-evolution models including the evolution of dust content have been 
developed (Wang 1991; Takagi & Arimoto 1997; Lisenfeld & Ferrara 1998, hereafter LF98; 
Dwek 1998, hereafter D98). The processes of dust formation and destruction by SNe are 
taken into account in LF98, in order to explain the relation between dust-to-gas mass ratio 
and metallicity of dwarf irregular galaxies and blue compact dwarf galaxies. We can find 
detailed mathematical frameworks to calculate the dust abundance in any galactic system 
in D98, in which the accretion of heavy elements onto preexisting dust grains is considered 
in addition to the processes considered in LF98. 

In this paper, we examine the dust-to-gas ratio of nearby star-forming galaxies with a 
simple model. We use spiral galaxies as the sample of star-forming galaxies. First of all, in 
the next section, we apply instantaneous recycling approximation to the model equations in 
D98 to describe a simplified relation between dust-to-gas ratio and metallicity. Our model 
is a one- zone model (i.e., the spatial distribution of gas, metal, and dust within a galaxy 
is not taken into account). The result is compared with observational data in the same 
section. Finally, we present summary and conclusion in §3. 



2. MODEL 

2.1. Basic Equations 

We, here, present the equations which describe the changing rate of the mass of gas, 
metal and dust in a galaxy. Here, we adopt a simple one-zone model, because we are 
interested in the global properties of galaxies. For the model treating the radial distribution 
of gas, element, and dust in a galaxy, see D98. After some modifications, equations (6)-(8) 
of LF98 are written as 

+ (1) 

-X^ + Ei, (2) 

M dl {\ - fj) M dl 
finiEi- afiXiip H . (3) 

7~acc 7"SN 

Here, M g is the mass of gas. The metal is labeled by % (i — O, C, Si, Mg, Fe, • • •), and Mi 
and M di denote the total mass of the metal i (in gas and dust phases) and the mass of 
the metal i in the dust phase, respectively. The star formation rate is denoted by E is 
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the total injection rate of mass from stars; Xj is the mass fraction of the element % (i.e., 
Xi = Mi/Mg); Ei is the total injection rate of element % from stars; fa is the mass fraction of 
the element i locked up in dust (i.e., fa = M di /Mi). The meanings of the other parameters 
in the above equations are as follows: fa n i represents the dust mass fraction in the injected 
material, in other words, the dust condensation efficiency in the ejectaQ; a refers to the 
efficiency of dust destruction during star formation (a = 1 corresponds to destruction of 
only the dust incorporated into the star, and a > 1 [a < 1] corresponds to a net destruction 
[formation] in the star formation); r acc is the accretion timescale of the element i onto 
preexisting dust grains in molecular clouds; t sn is the timescale of dust destruction by SN 
shocks. We hereafter assume that only dust incorporated into stars is destroyed when the 
stars form; i.e., a = 1. 

Our equations above are different from the equations in LF98 in the following two 
ways. The one is that we neglect the outflow term (i.e., we set W = in the equations of 
LF98), since we expect that outflows of gas from spiral galaxies, in which we are interested, 
are not effective because their gravitational potentials are much deeper than those of dwarf 
galaxies, which LF98 treated. For this reason and for simplicity, we here neglect the outflow 
terms. The other is the addition of the third term of the left-hand side in equation @. 
This term represents the accretion of element % onto the preexisting dust grains. The 
process is effective in molecular clouds, where the density of gas is high enough for the dust 
particles to grow through collisions with the atoms of heavy elements. The timescale r acc 
is determined from the lifetimes of molecular clouds. The detailed discussion about the 
accretion process can be found in §7 of D98. 

Here, we adopt the instantaneous recycling approximation (Tinsley 1980): Stars less 
massive than m\ live forever and the others die instantaneously. This approximation allows 
us to write E and respectively, as E = TZip and E^ = (TZXi + y^ijj, where TZ is the 
returned fraction of the mass that has formed stars which is subsequently ejected into the 
interstellar space, and 3^ is the mass fraction of the element i newly produced and ejected 
by stars.0 1Z and 3^; can be obtained using the following formulae (Maeder 1992): 



In equation ([|), <j>(m) is the initial mass function (IMF), and m u is the upper mass cutoff 

In this formalism, we assume that this efficiency in stellar winds is the same as that in SN ejecta. The 
different efficiencies in the various ejecta are treated in D98. 

2 TZ = R and y i = y(l - R) for the notations in LF98. 
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(4) 




(5) 
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of stellar mass. The IMF is normalized so that integral of m<j)(m) in the full range of stellar 
mass becomes 1. Thus, 4>(m) has dimension of the inverse square of mass. In equation (|5|), 
w m is the remnant mass (w m = 0.7M & for m < 4M Q and w m = 1.4M for m > 4M ) and 
Pi(m) is the mass fraction of a star of mass m converted into the element i and ejected. 

Using the above parameters TZ and y iy equations become 

1 dM„ 



if) dt 
M g dXi 

ip dt 
M g dVi 

if) dt 



-1+K, (6) 

MnXi + yi)-[a-i + n- p m {i - /*) + /3 sn ]a, (8) 



where = M di /M g = fcXi, and /5 acc and /5 S n are, respectively, defined by (3 acc = M g /(r acc if)) 
and Psn = Mg/^TsNip)- (See also Eqs. 10-12 in LF98.) We can regard /?sn as constant 
in time (§6.2 of D98). The reason is as follows. The destruction rate of dust by SNRs is 
written as (Mdj/M g )eMsNR-RsN, where e is the grain destruction efficiency, Msnr is the 
total mass of interstellar gas swept up by a SNR during its lifetime, and i?sN is the number 
of SNe (Type I and Type II) per unit time. Since this should be equivalent to the last 
term of right-hand side of equation (|3|), r SN = M g /(eM SNR i? SN ). Because the dependence of 
eM SNR on the gas density is weak (Jones, Tielens, & Hollenbach 1996), we regard eM SNR 
as a constant. If we reasonably assume that i?sN is proportional to the star formation rate, 
we derive tsn oc M g /if>. This is why /5sn = ^/(tsnVO can ^ e treated as a constant. The 
Galactic value shows f3$N ~ 5 (LF98). This value corresponds to r SN ~ 10 8 yr, which is 
consistent with Jones et al. (1994) and Jones, Tielens, & Hollenbach (1996). The relation 
race ^ r SN /2 (D98) leads to /3 acc ~ 2f3 SN ~ 10. 

Combining equations (^) and (§), we obtain the following differential equation of T>i as 
a function of Xf. 

dV 3 V 2 

y iT ^ = fUKXt + y i )-(a-l + K- 3 acc + f3 SN )T>i - (9) 

where we used the relation of fa = Vi/Xi. 



2.2. Comparison with Observations 

The result of the numerical integration of equation (0) is compared with observational 
data of nearby spiral galaxies. We take again i = O (see also §2.1). Moreover, we make 
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an assumption that total mass of dust including both oxygen dust and carbon dusif] is 
proportional to that of oxygen in the dust phase. In other words, 

V = YPi = GV , (10) 

i 

where C is assumed to be constant for all spiral galaxies. We adopt the Galactic value, 
C ~ 2.2 according to Table 2.2 of Whittet (1992). 

For further quantification, we need to fix the values of 71 and 3^ for the traced element 
(i — O). We choose oxygen as the tracer, because (i) Most of it is produced in type II SNe 
which are also responsible for the shock destruction of dust grains; (ii) oxygen is the main 
constituent of dust grains (see also LF98). The point (i) means that the instantaneous 
recycling approximation may be a reasonable approximation for the investigation of oxygen 
abundances, since the generation of oxygen is a massive-star- weighed phenomenon. In 
other words, results are insensitive to the value of m\ defined in equations (|) and (Q). 
According to LF98, we put m\ = 1M and m u = 12OM . We use a power-law form of the 
IMF: <f)(m) oc m~ x . Two kinds of the IMFs are investigated: x = 2.35 (Salpeter 1955) and 
x = 2.7 (Scalo 1986). If we assume that the lower cutoff of the stellar mass is equal to 
TOi(= 1M Q ), (K, y ) = (0.79, 1.8 x 10~ 2 ) and (0.68, 9.6 x 10~ 3 ) for the IMFs of Salpeter 
and Scalo, respectively (LF98). Adopting the lower cutoff of O.1M instead, the values 
become (K, y ) = (0.32, 7.2 x 10~ 3 ) and (0.13, 1.8 x 10~ 3 ) for the Salpeter and Scalo 
IMFs, respectively. We use the former values for 71 and 3fo- We treat the latter values only 
in 82.2.3. 



2.2.1. Dependence on /3 acc 

First, we examine the dependence of Xq- T> relation on the parameter /3 acc , the 
efficiency of accretion of heavy element onto the preexisting dust grains. The result is 
shown in Figure 1. 

The data points are from Issa, MacLaren, & Wolfendale (1990; see also the references 
therein). The dust-to-gas ratio is calculated from the interstellar extinction and column 
densities of H I and H2. Each value is evaluated at a galactocentric distance equivalent to 
that of the Sun (the galactocentric distance of each galaxy is properly scaled according to 
its size). 



3 We should note that the carbon dust is not traced by oxygen. 
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We show the cases of /3 acc = 5, 10, and 20 (the solid, dotted, and dashed lines, 
respectively), with the relation /3 acc = 2/3 SN fixed. The increase of /3 acc means that accretion 
of heavy elements onto preexisting dust becomes efficient. Thus, for a fixed value of the 
metallicity, dust-to-gas ratio increases as /3 acc increases. We also show the case of /5 acc = 
and /?sn = 10 (long-dashed line), which is far from reproducing the data. This means that 
we should take into account the accretion onto the preexisting dust as Dwek & Scalo (1980) 
asserted. Thus, the models in LF98, which do not include the term of the accretion, cannot 
be applied to spiral galaxies. The equations in LF98 can be applied to dwarf galaxies, where 
the accretion process may be inefficient because of their low metallicity. The accretion 
process is properly considered in D98. 



2.2.2. Dependence on f ini 

In Figure 2, the solid, dotted, and dashed lines show the V-X relations for / in0 = 0.4, 
0.1, and 0.01, respectively (in the case of /5 acc = 2/?sn = 10 and the Salpeter IMF). In LF98, 
/ino — 0.4 is suggested to be an optimistic upper limit. Figure 2 shows that the larger 
efficiency of metal production from heavy elements leads to the larger dust-to-gas ratio. 

In the limit of X — > 0, the solution of equation @ reduces to V ~ Cf in oX ( see 
also LF98), so that V scales linearly with / in0 . Thus, the precise value of / in o can be 
constrained by low-metal galaxies. Indeed, the difference among the three cases is significant 
in the low-Xo region in Figure 2. It seems unlikely that f\ n o 5s 0.01, since the data of the 
low-metal spirals cannot be reproduced by such small /mo 



2.2.3. Dependence on IMF 

The difference of mass function is represented by different 71 and 3^o (§2.1). We 
examine how the difference affects the result. We present the result in Figure 4. The 
Salpeter and Scalo IMFs are represented by the solid and dotted lines, respectively (lower 
ones). As commented above, if we use the lower mass cutoff for the stellar mass, 71 and 
3^o become smaller. We also show the case of lower mass cutoff (O.1M ; the upper line 
for each IMF). The differences among the results are too small (comparable to the typical 
observational errors) to determine which line is preferable. 
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3. SUMMARY AND CONCLUSIONS 

Based on the models proposed by LF98 and D98, we have examined the dust 
content in nearby spiral galaxies. We simplified the basic equations in D98 and applied the 
instantaneous recycling approximation to them. We treated one-zone model equations which 
describe the changing rate of dust-to-gas ratio including the terms of dust condensation 
from heavy elements ejected by stars, destruction by SN shocks, destruction in star-forming 
regions, and accretion of elements onto preexisting dust grains (§2.1). Assuming that the 
total dust mass is proportional to the mass of oxygen in the dust phase, we compared 
the results^ with the data of nearby spiral galaxies (for the detailed comparison with 
the Galactic values, see D98, in which the data are excellently reproduced). Finally, we 
summarize the conclusions reached in this paper: 

(i) Unless we do not take into account the accretion process of heavy elements onto 
the preexisting dust particles, we cannot explain the observed relations between dust-to-gas 
ratios and metallicities of nearby spiral galaxies. The accretion process is important for 
high-metallicity systems, since the collision rate between heavy-element atoms and dust 
grains is large. The data of dwarf galaxies, which are generally low-metallicity systems, 
may be explained even if the accretion process is neglected (LF98). The recent model by 
D98 properly treated the accretion. 

(ii) The efficiency of dust production from heavy element (denoted by /i n j) can be 
constrained by the galaxies with low metallicity. It is unlikely that f ini <^ 0.01. 

(iii) As for IMFs, both the Salpeter and Scalo types are consistent with the observed 
relation between dust-to-gas ratio and metallicity. 

We are grateful to S. Mineshige for continuous encouragement. We wish to thank E. 
Dwek for many comments which improved discussions in this paper and our knowledge of 
recent studies on dust. We acknowledge the Research Fellowship of the Japan Society for 
the Promotion of Science for Young Scientists. 
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FIGURE CAPTIONS 

FIG. 1 — The relations between dust-to-gas mass ratio (V) and oxygen abundance (oxygen 
mass fraction X Q ) for various /9 acc and Psn (dust accretion efficiency onto the preexisting 
dust grains and dust destruction efficiency in the SN shocks, respectively), with the relation 
f3 &cc = 2/9sn kept. The data points for nearby spiral galaxies are from Issa, MacLaren, & 
Wolfendale (1990). Here, / in0 is fixed to 0.1, and the Salpeter IMF is adopted. The solid, 
dotted, and dashed lines represent the cases of (3 acc = 2(3sn = 5, 10, and 20, respectively. 
The long-dashed line shows the case of no accretion process onto the preexisting dust grains 
(P^c = and (3 SN = 5). 

FIG. 2 — The same as Fig. 1 but for various / in0 (condensation efficiency of oxygen into 
dust grains). The other parameters are set to /3 acc = 2f3 SN = 10. The Salpeter IMF is 
adopted. The solid, dotted, and dashed lines represent the cases of / in o = 0.5, 0.1, and 
0.01, respectively. 

FIG. 3 — The same as Fig. 1 but for the different IMFs (Salpeter and Scalo; solid and 
dashed lines, respectively). Two lower cutoffs of stellar mass are investigated for each IMF: 
the cutoff of 1M Q and 0.1M© (the lower and upper lines, respectively). The parameters are 
fixed to /5 acc = 2/Ssn = 10 and / in0 = 0.1. 



